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PREFACE 


The  Amy  Symposium  on  |  So  lid  Mechanics,  1986,  represents  the 
tenth  in  the  series  of  bienniel  meetings  of  researchers  from 
government,  academia  and  industry  working  in  the  field  of  solid 
mechanics  in  support  of  Army,  Navy  and  Air  Force  systems. 

This  document  contains  the  extended  abstracts  of  work-in- 
progress  selected  for  presentation  at  the  1986  symposium  and 
inclusion  in  the  proceedings.  A  companion  document  contains  the 
full  technical  manuscripts  of  completed  work  selected  for  presen¬ 
tation  at  the  meeting. 

Initially  organized  in  1966  to  answer  the  needs  of  Army 
researchers  for  a  forum  for  presentation  and  critical  review  of 
work  being  accomplished  by  Army  laboratories  and  their  contrac¬ 
tors,  the  symposium  was  expanded  in  1968  to  include  the  partici¬ 
pation  of  the  Navy,  Air  Force  and  other  government  agencies  active 
in  areas  of  Amy  interest.  During  the  ensuing  years  the  symposium 
fulfilled  its  original  purpose  by  providing  an  excellent  opportu¬ 
nity  for  the  interaction  of  mechanics  researchers  in  the  three 
services  as  well  as  the  Departments  of  Energy  and  Transportation. 

For  each  meeting  the  symposium  committee  selects  a  central 
theme  representing  a  technical  or  system-related  area  of  current 
major  interest  to  the  Army.  The  theme  "Lightening  of  the  Force" 
selected  for  this,  the  tenth  symposium,  represents  a  major  Amy 
priority  in  order  to  enhance  both  fighting  capability  on  the 
battlefield  and  the  capability  for  rapid  deployment  of  fighting 
forces  to  remote  locations.  The  full  papers  and  extended  ab¬ 
stracts  selected  for  presentation  and/or  inclusion  in  the  proceed¬ 
ings  are  representative  of  the  extensive  research  and  advance  de¬ 
velopment  being  carried  on  in  support  of  this  central  theme. 

Although  centered  on  the  applied  research  (6.?)  area,  the 
area  of  interest  of  the  symposium  spans  the  range  from  basic  (6.1) 
research  through  advanced  development  (6.4).  The  papers  selected 
for  this  symposium  provide  this  broad  mixed  coverage  of  the  cen¬ 
tral  theme. 

On  this  tenth  anniversary  a  review  of  the  record  of  the  Army 
Symposium  on  Solid  Mechanics  shows  that  it  is  an  appropriate  and 
popular  event  that  has  filled  an  important  need  over  the  twenty 
years  of  its  existence.  It  is  fitting  now  to  acknowledge  the 
foresight  and  guidance  of  the  individuals  who  initiated  and  stead¬ 
fastly  supported  the  symposium  from  the  beginning.  Therefore  to 
Joseph  I.  Bluhm,  who  throughout  his  career  was  the  leading  advo¬ 
cate  of  solid  mechanics  research  and  engineering  in  the  Army;  and 
to  the  late  Dr.  John  M.  Slepetz  who  led  and  guided  the  Army 
Materials  Technology  Laboratory  and  the  Army  to  and  through  the 
intricacies  of  advanced  composites. 
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COMPUTED  NEUTRON  TOMOGRAPHY 
AND 

CODED  APERTURE  HOLOGRAPHY 
FROM  REAL  TIME  NEUTRON  IMAGES 


* 

MARK  F.  SULCOSKI 
University  of  Virginia 
P.  0.  Box  7432 

Charlottesville*  Virginia  22906 


EXTENDED  ABSTRACT 


Neutron  radiography  is  well  known  as  a  non-destructive  examination 
(NDE)  technique  for  low  atomic  number  materials.  The  increased  interest 
in  lightweight  composite  materials  using  low  atomic  number  elements 
(such  as  carbon)  require  modern  imaging  techniques  such  as  computed 
neutron  tomography  and  coded  aperture  neutron  holography  to  complement 
other  NDE  techniques.  The  unique  interaction  mechanisms  of  neutron 
radiation  as  compared  to  x-rays  allow  higher  detection  structures,  adhe¬ 
sive  bonding  layers  in  laminated  structures  and  delamination  detection 
in  composite  structures. 

The  purpose  of  the  continuing  work  described  herein  is  to  develop 
and  investigate  the  uses  of  neutron  tomography  and  holography  for  NDE 
applications.  Much  of  the  early  work  in  neutron  tomography  was  limited 
by  slow  and  tedious  methods  of  optical  densitometer  scanning  of  static 
neutron  radiographs  to  obtain  data.  The  use  of  boron  trifluoride  coun¬ 
ters  rather  than  film  improved  the  usefulness.  This  paper  describes 
the  use  of  a  real  time  neutron  imaging  system  coupled  with  a  dedicated 
image  processing  system  to  obtain  neutron  tomographs. 

Experiments  to  date  have  utilized  a  Thomson-CSF  neutron  camera 
coupled  to  a  Digital  Equipment  Corporation  LSI  11/23  based  system  used 
for  image  processing.  Experiments  have  included  the  following: 

Configuration  of  a  reactor  neutron  beam  port  for  neutron 
i mag i ng . 

Development  and  implementation  of  a  convolution  method 
tomographic  algorithm  suitable  for  neutron  imaging. 


Results  to  date  have  demonstrated  the  proof  of  principal  of 
this  neutron  tomography  system: 


Tomographs  from  real  time  images  of  1.6  cm  (0.63  in)  wooden 
test  objects  have  achieved  near  theoretical  resolution  of 
1.3  mm  <0.051  in)  voids  using  S00  projections  with  no  image 
processing . 

Computation  time  of  approximately  30  seconds  per  projection. 

Preliminary  benefits  of  image  processing  before  tomography 
data  acquisition  investigated. 


Future  work  is  planned  to  include  the  following: 

Investigation  of  alternate  algorithms  to  reduce  the  number 
of  projections  necessary  for  a  given  resolution. 

Investigation  of  combinations  of  image  processing  techniques 
and  other  data  preprocessing  algorithms  to  allow  larger 
objects  to  be  imaged  and  to  increase  resolution. 

More  studies  of  void  detection  and  delamination  detection 
for  various  materials  and  structures. 


Coded  aperture  neutron  holography  is  under  investigation  using  a 
cadmium  Fresnel  zone  plate  as  the  coded  aperture  and  the  real  time 
imaging  system  as  the  detection  and  holographic  reconstruction  system. 
Coded  aperture  imaging  utilizes  the  zone  place  to  encode  scattered 
neutron  radiation  from  the  entire  volume  of  an  object.  The  scattered 
radiation  pattern  recorded  at  the  detector  is  used  as  input  data  to  a 
convolution  algorithm  which  reconstructs  the  scattering  source  (object) 
Thus,  ideally,  all  the  tomography  planes  of  the  object  are  recorded 
simultaneously  allowing  a  three  dimensional  neutron  image  of  the  object 
This  technique  has  not  yet  been  successfully  implemented  and  is 
still  under  development.  Current  ongoing  work  includes: 

Investigation  of  techniques  to  increase  object  scattered 
neutron  flux  and/or  increased  detector  sensitivity. 

Comparison  of  several  methods  to  improve  signal-to-noise 
ratio  of  the  neutron  detection  system. 


# 

Present  address:  U.  S.  Army  Foreign  Science  and 

Technology  Center 
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ABSTRACT 


Resisting  strength  and  other  properties  of  engineering  materials  used 
for  the  production  of  laminated  orthotropic  composites  arc  known  to  have 
considerable  spatial  variations.  In  the  present  study,  a  particular  emphasis 
is  placed,  however,  on  the  probabilistic  model  for  spatial  variation  of  mate¬ 
rial  strength  and  on  the  corresponding  statistical  size  effect.  The  principal 
idea  lies  in  the  interpretation  that  the  material  strength  can  be  idealized  as 
a  multidimensional  random  field.  This  interpretation  together  with  the  re¬ 
cently  developed  method  of  digital  simulation  of  multidimensional  random 
helds  by  use  of  ARMA  models  makes  it  possible  to  demonstrate  the  statis¬ 
tical  size  effect  in  terms  of  numerical  examples.  Tha  ARMA  model  method 
mentioned  above  uses  a  recursive  equation  whose  coefficient  matrices  are 
determined  in  accordance  with  the  prescribed  autocorrelation  function  of 
the  multidimensional  random  field.  Using  the  recursive  equation  with  these 
coefficient  matrices,  we  can  generate  with  substantial  computational  ease, 
sample  functions  of  the  field  over  a  large  two-  or  three-  dimensional  do¬ 
main.  An  alternative  method  for  the  digital  simulation  of  multidimensional 
random  fields  is  th°  well-known  Fast  Fourier  Transform  (FFT)  method.  It 
should  be  pointed  out  however,  that  the  ARMA  model  method  is  consider- 
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ably  more  efficient  in  terms  of  computational  time. 

Consider  now  a  specific  example  of  a  laminated  orthotropic  compos¬ 
ite  plate  consisting  of  two  laminae.  Each  lamina  is  divided  into  an  equal 
number  of  elements  (4  •  6  =  24  for  the  case  considered  -  see  Fig.  1). 


Fig.  1  Division  of  plate  into  48  elements 

Therefore  the  whole  plate  is  divided  into  2  •  24  =  48  elements.  We 
can  now  assign  values  for  the  material  strengths  in  the  three  principal  di¬ 
rections  (longitudinal,  transverse  and  shear  strengths  respectively)  for  each 
one  of  the  48  elements.  This  can  be  achieved  by  considering  that  for  each 
lamina  the  material  strengths  in  the  three  principal  directions  consist  a 
three-variate  two-dimensional  (3V-2D)  random  field.  In  this  way,  each  one 
of  the  48  elements  is  assigned  with  a  different  triad  of  values  for  the  material 
strengths  in  the  three  principal  directions. 

Considering  now  uniaxial  stress  of  the  laminated  orthotropic  composite 
plate,  we  can  start  from  a  very  small  stress-level  a  and  check  whether  one  or 
more  of  the  48  elements  fail  under  the  given  stress  a.  The  Tsai-Hill  failure 
criterion  shown  below  in  Eq.  (a),  can  be  used  to  perform  this  check.  It 
should  be  noted  that  this  failure  check  is  performed  48  times,  once  for  each 
element.  The  Tsai-Hill  failure  criterion  is  given  by  : 
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X2  VS2  x2’  Y2  a2  v/ 

where  X,  Y,  S  are  the  strengths  in  the  three  principal  directions  (longitu¬ 
dinal,  transverse  and  shear  respectively),  a  is  the  externally  applied  axial 
stress  and  0  is  the  angle  between  the  principal  direction  X  and  the  direction 
of  application  of  the  uniaxial  stress  a.  Of  course,  any  other  failure  criterion 
can  be  used  instead  of  the  Tsai-Hill  one,  by  making  slight  modifications 
to  the  above  procedure.  For  example  some  other  failure  criterion  might 
involve  only  two  independent  strengths  instead  of  the  three  of  the  Tsai-Hill 
criterion.  In  this  case,  a  two-variate  two-dimensional  random  field  has  to 
be  generated. 

If  no  failure  is  detected  at  the  original  stress-level,  the  stress-level  can 
be  increased  by  steps  up  to  the  point  where  some  element(s)  fails  for  the  first 
time.  The  value  of  the  externally  applied  uniaxial  stress  at  this  level  can 
be  considered  as  the  strength  of  the  whole  laminated  orthotropic  composite 
plate  in  uniaxial  tension. 

Using  Monte  Carlo  techniques  the  above-described  procedure  can  be 
repeated  a  number  of  times  for  different  realizations  of  the  3V-2D  random 
field  describing  the  strengths  of  the  laminated  composite  plate  in  the  three 
pricipai  directions.  In  this  fashion  several  different  values  for  the  strength  of 
the  whole  laminated  plate  can  be  obtained.  Then,  the  distribution  function 
of  those  values  can  be  estimated  by  the  use  of  some  test  of  goodness  of  fit 
(e.g.  probability  paper). 

Finally,  in  order  to  examine  the  statistical  size  effect,  the  following 
procedure  can  be  followed.  The  same  problem  as  the  one  described  abr*ve 
to  find  the  strength  of  the  whole  laminated  composite  plate  is  solved  again 
by  considering  in  the  beginning  twice  the  area  of  the  original  laminated 
plate,  then  thrice  and  so  on,  however  always  using  the  same  autocorrela¬ 
tion  function  for  the  3V-2D  random  field  describing  the  material  strength 
variability  in  the  three  principal  directions.  In  this  way,  conclusions  can  be 
drawn  about  the  statistical  size  effect  due  to  the  spatial  strength  variation. 


A  DIRECT  METHOD  TO  OBTAIN  CONSTITUTIVE  MODEL  CONSTANTS 
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GORDON  R.  JOHNSON 
Principal  Engineering  Fellow 
Honeywell  Defense  Systems  Division 
Edina,  MN  55436 


INTRODUCTION 


The  cylinder  impact  test  was  introduced  by  G.  I.  Taylor  in  1948.  Since  that  time 
the  merits  of  this  test  have  been  debated.  Some  would  argue  that  this  test  should  not 
be  used  to  obtain  constitutive  data  because  a  wide  range  of  unknown  strains,  strain 
rates,  temperatures,  and  stresses  vary  throughout  the  test  specimen,  and  throughout 
the  duration  of  the  test.  It  is  better,  they  would  say,  to  obtain  data  from  dynamic  tension, 
compression,  and/or  torsion  tests,  in  which  the  variables  are  better  defined,  and  then 
use  the  cylinder  impact  test  to  provide  an  independent  check  of  the  data. 

Others  would  argue  that  dynamic  tension,  compression,  and/or  torsion  tests 
require  expensive  test  facilities  and  that  these  tests  are  limited  in  their  ability  to  achieve 
the  range  of  strains  and  strain  rates  of  interest.  Conversely,  the  cylinder  impact  test 
is  simple,  inexpensive,  and  provides  data  for  the  appropriate  range  of  strains  and  strain 
rates.  During  the  past  twenty  years,  several  researchers  have  used  this  test  to  esti¬ 
mate  dynamic  flow  stresses  for  various  materials. 

The  work  reported  in  this  paper  is  directed  at  developing  an  explicit  method  (as 
opposed  to  previous  trial-and-error  methods)  to  obtain  constants  for  computational  con¬ 
stitutive  models. 

DISCUSSION  AND  RESULTS 

Figure  1  shows  a  comparison  of  cylinder  impact  test  results  and  corresponding 
EPIC-2  computational  results  for  OFHC  copper.  The  test  data  are  from  Reference  1 
(Johnson  and  Cook,  Proceedings  of  the  Seventh  International  Symposium  on  Ballistics, 
1983).  The  adiabatic  stress-strain  relationships  are  given  in  Figure  2  and  the  strength 
model  constants  used  for  the  four  cases  are  given  in  Table  I.  The  adiabatic  stresses 
are  shown  only  to  the  maximum  strains  attained  in  the  computed  results. 
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NOTES: 


•  Of  HC  COPPER  CYLINDER  (L„=  26.4mm.  D„=  7.6mm,  L  =  16.2mm,  D  a  13.6mm) 


•  IMPACT  VELOCITY  *  190  m/« 

•  TEST  RESULTS  INDICATED  BY  DOTS •  ••  • 

•  EQUIVALENT  PLASTIC  STRAIN  CONTOURS  SHOWN  AT  INTERVALS  OF  0.6 


Figure  1.  Comparison  of  Test  Results  with  Computed  Results  Using  Various 
Strength  Model  Constants 


CASE  B 


-CASE  D  (103<  C  <  105> 


CASE  C 


CASE  A 


EQUIVALENT  PLASTIC  STRAIN,  C 


Figure  2.  Adiabatic  Stress-Strain  Relationships  for  the  Various  Strength  Models 


Table  I.  Strength  Model  Constants  for  the  Computed  Results 
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98 

469 

368 

.66 
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0. 

.026 

1.00 
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For  Case  A  the  length,  L,  of  the  deformed  cylinder  is  matched  with  the  computa¬ 
tional  result.  Because  only  one  deformed  dimension  is  matched,  only  one  independent 
strength  constant  (representing  a  constant  flow  stress)  can  be  obtained.  It  can  be  seen 
that  there  are  significant  discrepancies  between  the  test  and  computational  results 
at  the  deformed  end  of  the  cylinder. 

Case  B  provides  for  linear  strain  hardening.  By  matching  both  the  deformed 
length,  L,  and  the  maximum  diameter,  D,  it  is  possible  to  obtain  the  two  constants 
for  the  linear  hardening.  Here  the  computed  result  is  in  good  general  agreement  with 
the  test  result.  If  this  model  would  be  applied  to  larger  strains  than  those  experienced 
in  the  computed  results  of  Case  B,  however,  it  would  probably  lead  to  excessively  high 
stresses.  It  is  well  known  that  thermal  softening  tends  to  decrease  the  rate  of  strain 
hardening  at  large  strains.  In  general,  if  the  basic  model  accurately  represents  material 
behavior,  then  the  results  can  probably  be  extrapolated  beyond  the  data  base. 

Cases  C  and  D  are  based  on  the  following  model  for  the  von  Mises  tensile  flow 
stress  (Reference  1): 

o  =  [A  +  Be"][1  +  Cln  i*)[l  -  T*m]  (1) 

where  i  is  the  equivalent  plastic  strain,  i*  =  Hi0  is  the  dimensionless  plastic  strain 
rate  for  £0  =  1.0s-1,  and  T*  is  the  homologous  temperature.  The  five  material  con¬ 
stants  are  A,  B,  n,  C,  m.  The  first  three  constants  (A,  B,  n)  represent  the  strain  harden¬ 
ing  characteristics  of  the  material  and  are  the  ones  which  are  determined  from  the 
test  results.  Because  there  are  now  three  constants,  it  is  possible  to  match  three  deformed 
dimensions— the  length,  L,  the  maximum  diameter,  D,  and  an  intermediate  diameter, 
W.  The  strain  rate  constant,  C,  and  the  thermal  softening  constant,  m,  must  be 
approximated  or  obtained  from  other  sources. 


Case  C  is  the  result  obtained  if  nothing  is  known  about  the  strain  rate  or  ther¬ 
mal  softening  characteristics  of  the  material.  Here,  the  strain  rate  constant  is  set  to 
C  =  0,  and  the  thermal  softening  is  assumed  to  be  linear,  or  m  =  1.0.  Case  D  uses  strain 
rate  and  thermal  softening  data  from  Reference  1.  Even  though  there  is  little  appar¬ 
ent  difference  between  Cases  C  and  D  for  this  problem,  for  some  other  problems  (with 
wider  ranges  of  strain  rates  and  temperatures)  the  Case  D  constants  will  probably  pro¬ 
vide  better  results. 

METHOD  TO  OBTAIN  CONSTANTS 


The  method  consists  of  applying  the  following  equation  in  an  iterative  manner 
until  a  satisfactory  set  of  constants  is  attained: 


{Ax}  =  [K]{  Ay } 

| 

fAA> 

}  1 

'a.-) 

where  { Ax }  =  < 

\  AB  ( 

[  and  { Ay }  =  < 

|ad 

1 

f  An  1 

^  j 

1  i 

!  AW 

l  ) 

For  a  given  set  of  material  constants  (A,  B,  n)  it  is  possible  to  compute  the  deformed 
geometrical  dimensions  (L,  D,  W).  These  deformed  dimensions  are  compared  to  the 
dimensions  of  the  test  specimen  and  the  differences  are  expressed  as  AL,  AD,  and  AW. 
These  differences  can  then  be  used  to  find  the  corresponding  material  constant  changes 
which  will  lead  to  improved  material  constants.  This  is  accomplished  through  a 
3x3  operator  matrix,  [K],  where  Ky  is  the  change  in  material  constant,  Ax*,  due  to 
a  unit  change  in  the  deformed  geometry,  Ayj. 

The  operator  matrix,  [K],  can  most  easily  be  attained  by  first  developing  the  in¬ 
verse  to  this  matrix.  This  is  accomplished  by  running  a  series  of  four  computer  simu¬ 
lations  with  different  material  model  constants.  The  first  simulation  is  a  baseline 
computation  which  uses  an  initial  estimate  of  the  constants.  The  next  three  simula¬ 
tions  vary  each  of  the  three  constants  independently  (A',  B',  n').  This  procedure,  and 
the  generation  of  the  [K]  \  is  shown  in  Figure  3. 

The  operator  matrix  assumes  a  linear  relationship  between  incremental  material 
constant  changes  and  incremental  deformed  dimension  changes.  This  is  not  a  valid 
assessment  because  the  operator  matrix  is  actually  nonlinear.  However,  the  linear 
assumption  can  produce  accurate,  converging  results  if  the  four  simulations  used  to 
generate  the  operator  matrix  are  sufficiently  close  to  the  final  solution. 

For  the  limited  number  of  solutions  obtained  to  date,  convergence  was  experienced 
after  several  iterations.  In  these  cases,  the  initial  operator  matrix  was  used  for  all  iter¬ 
ations.  The  rate  of  convergence  can  eventually  be  increased  by  implementing  improved 
techniques  that  continually  update  the  operator  matrix  after  every  simulation. 
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SIMULATION 

MATERIAL 
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Figure  3.  Summary  of  Computer  Simulations  Required  to  Generate  Inverse  of 
Operator  Matrix 


SUMMARY  AND  CONCLUSIONS 

It  has  been  demonstrated  that  it  is  possible  to  directly  obtain  constitutive  model 
constants  from  cylinder  impact  test  results.  Future  work  will  be  directed  at  improving 
the  convergence  of  the  method,  and  applying  the  method  to  other  models  and  materials. 
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ABSTRACT 

Met  a 1 1 oqraphi c  studies^  have  shown  that  crack  propagation 

occurs  by  a  plastic  shear  instability  mechanism  in  high  strength 

martensitic  alloy  steels.  On  the  microscale,  propagation  occurs 

in  a  non-planar  fashion  following  bands  which  develop  along  the 

shear  directions.  Voids  evidently  play  a  major  role  in  this 

fracture  orocess,  with  microcracks  developina  between  interacting 

2 

voids.  In  shear  tests  of  high  strength  A3A0  steel  ,  voids  have 
been  found  within  shear  bands  at  debonded  sub-micron  sized  qrain 
refinement  particles.  It  has  been  speculated  that  the  microcracks 
linking  such  voids  form  as  a  result  of  the  coalescence  of  voids 
which  nucleate  at  smaller  scale  strengthening  particles.  It  is 
likely  that  this  secondary  void  nucleation  results  from  the  strain 
intensification  caused  by  the  dominant  voids.  Fracture  is 
expected  after  a  critical  number  of  microcracks  are  forme1 

The  present  mechanics  study  is  being  conducted  in  suooor 
an  alloy  development  research  program  attempting  to  oroducf 
ultra-high  strength  steel  alloy  (300-35C  ksi)  with  fracture 
touqhness  sufficiently  high  (6P-12G  ksi\/"in)  to  allow  a  wide  range 
of  structural  applications.  Success  of  this  program  will  provide 
an  important  new  material  selection  option  for  lightweight 
fracture  resistant  metallic  structures.  The  mechanics  work  is 
directed  to  establishing  a  quantitative  understanding  of  the 
void-softening'  mechanism  in  shear.  This  information  will  be  used 
to  guide  alloy  design  to  achieve  low  shear  bandinq  propensity. 

The  key  to  understanding  how  microscopic  voids  deteriorate  a 
plastically  deforminq  metal  is  knowledge  of  the  void  reqion 
stress/  strain  elevation  over  the  uniform  levels  that  would  exist 
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in  the  absence  of  voids.  Particularly  useful  would  be  knowledne 
of  the  void  region  state  when  the  remote  uniform  field  simulates 
the  deformation  occuring  in  a  band  of  localized  flow.  This  is  the 
analytical  problem  addressed  here.  A  set  of  elastic-plastic 
finite  element  analyses  has  been  conducted  to  establish  the  local 
field  between  a  pair  of  voids  close  enouqh  to  Droduce  significant, 
interacti ons . 

The  analyses  were  conducted  using  large  scale  comoutational 
facilities  and  an  elastic-plastic  finite  element  formulation  for 
non-hardening  P r and t  l -Reus s  constitutive  theory^.  The  formulation 
is  an  incremental  secant  stiffness  method  that  quarantees 
satisfaction  of  the  yield  condition  at  each  grid  point.  An 
adaptive  incrementation  scheme  is  employed  to  automate  load  step 
size  selection  and  to  control  load  path  discretization  error.  The 
formulation  traces  the  spread  of  elasticity  by  adjusting  element 
stiffnesses  as  the  imposed  displacement  is  incrementally 
increased.  The  mesh  used  to  model  the  void  pair  in  a  plane  strain 
deformation  band  is  shown  in  Figure  1.  There  are  approximately 
1 50C  nodes  (3000  degrees  of  freedom)  in  the  mesh  drawn.  However, 
symmetry  allows  half  of  the  model  in  simple  shear  analysis  and  a 
quadrant  for  uniaxial  extension.  All  four  sides  of  the  band  were 
constrained  according  to  nominal  deformation  mode.  In  simple 
shear,  the  edges  were  constrained  against  vertical  displacement, 
while  in  uniaxial  extension,  nodes  on  all  four  edges  were 
constrained  against  horizontal  motion.  Thus,  in  the  absence  of 
the  voids,  the  solutions  would  be  uniform  shear  (no  normal  strain) 
and  constant  uniaxial  extension  (a  single  non-zero  strain 
component ) . 

The  plastic  zone  development  for  the  simple  shear  case  is 
illustrated  in  Figure  2.  In  the  early  stages  of  loadinq,  plastic 
deformation  is  restricted  to  regions  at  the  void  surfaces,  rouqhly 
45  degrees  from  the  ligament.  As  these  surface  zones  develop,  a 
separate  distinct  plastic  zone  develops  along  the  ligament.  Once 
the  zones  link,  plastic  strain  localization  occurs  between  the 
voids.  Fiqure  3  illustrates  the  fact  that  over  the  early  stages 
of  loading,  straining  is  fairly  uniform  over  the  interior  of  the 
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ligament,  but  a  maximum  is  seen  to  develop  at  a  distance  of  P . 3  0 
ahead  of  the  voids.  The  value  of  the  maximum  local  shear  strain 
is  plotted  versus  nominal  strain  in  Figure  4.  The  rate  of 
increase  of  local  strain  with  respect  to  applied  strain 
continually  rises  with  the  final  intensification  rate  in  excess  of 
thirty  . 

The  plasticity  and  interaction  for  the  case  of  uniaxial 
extension  is  much  different  than  the  above  simple  shear  case.  The 
plastic  zone  at  a  uniaxial  strain  value  1.27  Y/F  is  shown  in 
Figure  5.  Also  shown  in  Figure  5  is  that  portion  of  the  plastic 
zone  that  satisfies  the  near  isocloric  condition  that  the 
magnitude  of  the  normal  strain  increments  differ  by  less  than  6?-'. 
Plasticity  spreads  from  the  void  surfaces  and  excessive  strain 
intensification  occurs  once  the  zones  merge.  The  numerical  data 
agrees  with  the  logarithmic  spiral  distribution  expected  from 
slipline  theory,  Figure  7.  Consistent  with  the  plastic  zone 
development,  the  maximum  strain  in  this  loading  case  is  seen  to 
occur  at  the  void  surfaces.  The  very  significant  intensification 
of  strain  suggests  that  substantial  void  growth  will  occur,  in 
agreement  with  known  results  for  the  growth  of  isolated  voids  in 
triaxial  stress  fields  . 

In  summary,  the  plasticity  results  for  the  test  cases  studied 
are  consistent  with  experimental  observations  of  void  sheet 
development  in  shear  and  void  growth  in  tension.  Further  work  is 
being  oursued  to  establish  interaction  features  over  a  ranne  of 
particle  distribution  parameters.  The  results  suaqest  that  the 
goal  of  establishing  alloy  development  guidelines  to  minimize 
shear  banding  propensity  using  micromechanics  is  achievable  and 
likely  to  be  very  worthwhile. 
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FIGURE  2  -  Development  of  plastic  instability  between  void 

pair  under  simple  shear  Plastic  zones  at  nominal 
shear  strain  levels  60.  77,  and  91*  of 
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A  Study  of  the  Crack  Tip  Parameters  that  Govern  Environmentally 
Assisted  Cracking  by  Use  of  the  Rising  Displacement  Test 


by 

R.  Mayville,  T.  Warren  and  P.  Hilton 
A. D .  Little,  Inc.,  Cambridge,  MA 


Abstract 

There  is  an  increasing  interest  in  the  use  of  rising  displacement  or  load 
tests,  similar  to  K-R  or  J-R  resistance  curve  tests,  to  characterize 
environmentally  assisted  cracking  (EAC)  in  metal/environment  systems. 

These  tests  have  been  used  to  determine  K  (K  )  and  to  study  the 

effect  that  crack  tip  parameters,  such  as  ,  dJ/dt  and  d5/dt,  have  on  EAC. 

This  investigation  is  an  experimental  study  of  the  'effect  a  wide  range  of 
applied  displacement  rates  has  on  EAC  for  AISI  4340  with  a q  -  1095  MPa  (159 
ksi)  in  flowing  synthetic  seawater.  Cantilever  beam  -  constant  load  - 
tests  were  performed  to  obtain  -  98  MPa  /m.  Rising  displacement 

tests  were  performed  over  four  orders  of  magnitude  with  25mm  thick  compact 
tension  specimens  under  ambient  conditions  and  in  environment  to  generate 
J-R  curves  and  to  determine  crack  velocities.  Under  ambient  conditions  the 
4340  exhibits  substantial  stable  tearing.  In  environment,  the  J-R  curves 
were  identical  to  the  ambient  J-R  curve  except  for  the  slowest  test,  which 
lasted  over  a  week.  The  J-R  curve  for  this  specimen  was  substantially 
lower  than  the  ambient  curve  and  KTT  converted  from  J  r  was  found  to  agree 
very  well  with  KIEAC-  The  curves  are  shown  in  Figure  1.  Two  experi¬ 
ments  were  also  performed  in  environment  in  which  the  displacement  -  stroke 
-  was  held  constant  after  some  stable  crack  extension  had  occurred.  The 
crack  velocities  from  these  and  the  other  compact  tension  tests  are  plotted 
in  Figure  2  against  K^.  and  indicate  that  the  a-Kj  relationship  is  unique 
for  environmental  cracking  for  (from  ambient  J^) .  In  addition, 

the  environmental  a-K^  curve  is 
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parallel  to  the  ambient  a-K^  curve  suggesting  that  the  ambient  mechanical 
cracking  mechanisms  play  a  major  role  in  environmental  cracking  for  this 
system.  This  is  supported  by  very  similar  ductile  rupture  fracture 
surfaces  for  the  two  cases.  This  work  has  been  supported  by  the  National 
Science  Foundation  and  by  the  Arthur  D.  Little  Corporate  Fund. 


Key  words  for  Army  Symposium:  fracture  mechanics,  experimental  methods, 

life  prediction. 
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Figure  2:  Crack  Velocity  Data  from  the  Compact  Tension  Specimens 
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EXPLICIT  AND  IMPLICIT  SOIL  MODELLING  FOR  DYNAMIC  SOIL-STRUCTURE 

INTERACTION  PROBLEMS 

Aaron  Das  Gupta,  Ph.D. 

US  Army  Ballistic  Research  Laboratory 
Aberdeen  Proving  Ground,  Maryland  21005-5066 


ABSTRACT 

Increased  application  of  nonlinear  viscoelastic  materials  in  structures 
and  consideration  of  operating  environments  require  modeling  of  interaction 
of  such  structures  with  the  surrounding  soil  medium  which  significantly 
influences  the  dynamic  response  when  subjected  to  an  impulse  or  blast  over¬ 
pressure  loading.  The  need  is  currently  accentuated  by  the  rapid  advances 
in  numerical  modeling  of  dynamic  problems  using  either  finite-difference 
or  finite  element  techniques  of  structural  analysis. 

While  various  models  are  available  to  describe  the  deformation  and  failure 
processes,  scarcity  of  dynamic  experimental  data  at  high  strain  rates  is  a 
serious  handicap  for  a  majority  of  representative  soil  materials.  As  a 
result,  many  computations  are  still  made  using  static  or  low  strain  rate 
data.  In  order  to  overcome  this  critical  shortcoming  in  soil-structure  inter¬ 
action  problems,  an  explicit  modeling  technique  could  be  improvised  by 
representing  the  soil  with  an  equivalent  mechanical  model^  consisting  of  a 
spring  and  a  dashpot  in  parallel  and  connecting  the  model  to  the  structural 
model  using  node-tie  element.  In  this  spring  support  condition,  the  base 
of  a  structure  in  contact  with  soil  could  be  flexibly  supported  using  an 
appropriate  value  of  the  spring  stiffness  while  the  top  of  the  structure 
could  be  loaded  directly. 

The  implicit  modeling  technique  where  nodal  tie  elements  are  used  to 
model  the  base  support  as  nonlinear  springs  for  shallow  buried  structures 
to  simulate  the  soil  medium  suffers  from  a  serious  shortcoming  in  sofaras 
the  indirect  loading  of  the  structure  through  shock  waves  passing  through 
the  soil  and  overpressure  loading  attenuation  effects  are  not  adequately 
modeled  resulting  in  unrealistic  soil-structure  interaction  and  erroneous 
results.  Hence  an  explicit  modeling  technique  is  recommended  whereby  a 
few  layers  of  soil  surrounding  the  base  of  the  structure  are  modeled  as 
an  assembly  of  elements  in  contact  with  the  structure.  The  interfacing 
cantact  elements  for  both  structure  and  soil  are  allowed  to  stick,  slide 
or  seperate  from  each  other  to  simulate  the  appropriate  boundary  condi¬ 
tions. 

The  explicit  modeling  technique  using  a  finite  element  code  require 
characterization  of  constitutive  relationship  of  geological  materials  from 
shock  compression  behavior  of  such  materials  and  their  Hugoniot  charac¬ 
teristics.  Although  the  equations-of-state  in  the  form  of  pressure- 
density  relationships  for  shock  propagation  in  such  mediums  can  be  non¬ 
linear,  the  characterization  technique  in  this  investigation  is  currently 
restricted  to  materials  with  linear  Hugoniot  characteristics  at  all 
pressure  ranges. 
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Three  representative  ^oil  mediums  i.e.,  saturated  and  unsaturated 
tuff,  oil  shale  and  sand  are  selected  for  this  study.  Under  the  assumption 
of  unidimensional  shock  propagation  and  conservation  law  of  mass,  momentum 
and  energy,  the  equations-of-state  could  be  transformed  to  bulk  moduli 
versus  volumetric  strain  relationship  using  the  methodology  described  in 
this  paper.  The  transformed  characterization  is  compatible  with  the 
tension  cutoff  curve  description  material  model  available  in  the  ADINA 
nonlinear  dynamic  finite  element  analysis  program.  However,  the  charac¬ 
terization  of  sand  is  valid  in  this  study  only  at  very  high  pressure 
regimes  where  sand  behaves  like  quartz  due  to  solid  phase  transition 
and  a  linear  Hugoniot  curve  fit  to  the  shock  velocity  versus  particle 
velocity  data  is  feasible. 

The  method  could  be  readily  extended  to  partially  saturated  soil  with 
nonlinear  characteristics  through  multilinear  approximation  during  the 
loading  and  unloading  phase.  The  resulting  constitutive  relationship  will 
facilitate  explicit  modeling  of  the  surrounding  medium  around  shallow- 
buried  structures  and  realistic  simulation  of  dynamic  soil-structurp 
interaction  which  can  significantly  alter  the  dynamic  response  behavior 
of  such  structures. 
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EXTENDED  ABSTRACT 

The  Mechanics  and  Structures  Branch  of  the  Interior  Ballistics  Division  of  the 
Ballistic  Research  Laboratory  has  embarked  upon  a  mission  program  tailored  to 
understanding  ^omposite  mater ia Island  their  potential  application  to  Army 
systems.  The  intent  of  the  program  is  to  provide  the  Amy  with  design  techniques 
and  prediction  methodologies  specifically  aimed  at  the  problems  peculiar  to  Army 
systems.  These  problems  include  high  rates  of  loading,  including  ballistic 
loadings  and  blast  loadings;  high  and  low  cycle  fatigue;  corrosion  and  erosion 
resistance;  and  the  necessity  to  predict  structural  response  of  thick  sections  to 
all  loading  environments.  The  intent  is  also  to  attempt  to  apply  new  and 
emerging  materials  developed  at  other  Amy  laboratories,  under  other  Amy 
programs,  and  those  materials  available  in  the  commercial  sector.  The  approach 
to  this  broad  scoped  and  general  program  has  been  to  decide  which  issues  are  the 
most  important  to  the  designers  of  some  specific  systems  and  components,  and 
initially  concentrate  upon  those  aspects  of  the  material  systems  which  are  not 
yet  well  understood.  Three  of  these  aspects  of  composite  material  behavior  have 
been  selected  as  the  most  critical  current  predictive  techniques  for  composite 
materials.  They  are  :  structural  response  to  high  rates  of  loading,  predicting 
structural  response  of  thick  sections  of  laid-up  laminated  structures,  and 
failure  criteria  for  thick  laminates.  The  remainder  of  this  abstract  will  focus 
on  the  work  in  progress  in  these  three  critical  areas. 


HIGH  HATE  LOADING 


We  are  developing  an  instrumented  projectile  that  can  be  launched  from  a  37mm 
gun  housed  in  an  indoor  range  at  BRL.  This  gun/projectile  system  can  be  used 
to  ballistically  test  a  sample  encased  within  the  instrumented  projectile.  By 
changing  the  charge  and/or  mass  and/or  geometry  of  the  projectile,  a  wide 
range  of  ballistic  loading  rates  and  states  of  stress  of  interest  can  be 
investigated.  It  is  intended  that  samples  of  composite  materials  will  be  made 
up  and  loaded  in  the  instrumented  projectile  for  tests.  Data  such  as  strain 
and  strain  rate  in  several  directions  and  at  several  locations  within  the 
sample  will  be  taken.  These  data  will  then  be  correlated  with  finite  element 
structural  simulations  of  the  tests  to  determine  the  efficacy  of  a  particular 
modelling  method,  and  also  to  benchmark  that  method.  In  this  way,  predictive 
techniques  of  structural  response  at  the  rates  of  loading  of  interest  can  be 
developed  and  benchmarked  before  being  applied  to  full  scale  Army  systems  and 
components . 

THICK  SECTIONS 

Army  structural  components  tend,  in  general,  to  be  larger,  thicker,  and  more 
massive  than  do  structural  components  for  agencies  such  as  the  Air  Force  or 
NASA.  The  major  studies  to  date  that  have  been  performed,  and  a  majority  of 
the  structures  which  have  manufactured  have  been  for  those  two  agencies. 

These  structures  tend  to  be  two  dimensional  curved  or  cylindrical  structures 
with  or  without  stiffening  ribs.  Army  Structures,  such  as  howitzer  components 
and  armored  vehicle  hulls,  tend  to  be  made  of  thick,  strong,  massive 
components.  To  date,  very  little  analytical  work  has  been,  performed  for  this 
type  of  structure. 

A  major  portion  of  the  composite  material  application  program  has  had  to  deal 
with  just  this  problem.  How  does  one  predict  structural  response  of  a  thick 
laminated  section.  Laminated  plate  theory  fairly  accurately  predicts  the 
response  of  thin  structural  components  such  as  wing  skins  or  rocket  motor 
housings,  but  the  theory  fails  to  accurately  predict  the  response  of  thick 
laminates.  We  have  therefore  taken  a  broad  approach  to  attacking  this 
problem.  First,  several  experts  from  within  and  outside  of  the  Army  labora¬ 
tories  have  been  contacted.  Experts  from  two  of  the  pre-eminent  centers  for 
composites  have  been  brought  under  contract  to  derive  failure  criteria.  This 
will  be  discussed  later  in  this  abstract. 

Secondly,  a  small  laboratory  is  being  set  up  which  is  specifically  geared 
toward  constructing  samples  of  thick  sections  for  testing.  Both  an  injection 
molding  machine  and  a  small  autoclave  will  be  used  to  make  up  test  samples  of 
thick  sections  for  mechanical  testing  and  benchmarking  of  structural  pre¬ 
dictive  techniques. 
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FAILURE  CRITERIA 


All  of  the  structural  response  predictions  that  can  be  made  are  worthless 
unless  we  can  predict  whether  or  not  the  structure  survives  its  loading 
environment.  One  of  the  most  important  parts  of  this  program  is,  therefore, 
the  prediction  of  failure,  and  the  criteria  whereby  that  failure  is  predicted. 
There  are  many  schools  of  thought  concerning  failure  prediction  in  laminated 
structures.  There  are  almost  as  many  as  there  are  experts  predicting  the 
failure.  They  can,  however,  be  broken  into  two  general  categories;  maximum 
stress  criteria  and  maximum  strain  criteria.  Experts  from  both  the  Air  Force 
Materials  Lab  and  Lawrence  Livermore  Laboratories  are  working  under  contract 
on  one  or  the  other  of  these  general  criteria.  It  is  intended  that  the 
results  of  these  two  studies  will  be  integrated  into  a  design  tool  for  use  by 
Army  systems  designers  which  will  allow  these  designers  to  accurately  predict 
the  response  and  failure  of  structural  components.  Enough  intelligence  will  be 
built  into  this  predictive  tool  to  allow  the  designer  to  select  from  several 
approaches  to  laying  up  a  laminate,  and  several  approaches  to  predicting  the 
response  and/or  failure  of  the  component  in  question. 


CLOSURE 

We  have  embarked  on  a  long  term  and  intentionally  aggressive  program  for  the 
development  of  predictive  techniques  and  design  methodologies  for  application 
of  composite  materials  to  Army  structural  components  and  systems.  The  program 
has  elements  of  pure  research  into  predictive  techniques  for  thick  sections 
and  failure  of  components,  and  elements  of  experimental  wprk  in  high  rates  of 
loading  common  to  Army  Ballistic  coraonents.  The  program  has  a  very  broad 
scope,  but  is  focused  upon  the  peculiar  needs  of  Army  component  and  systems 
designers.  It  meshes  well  with  ongoing  programs  at  other  major  Army  labora¬ 
tories  in  composite  materials  development,  and  it  also  has  brought  together  a 
large  body  of  knowledge  from  outside  experts  in  an  effort  to  capitalize  on  as 
much  of  the  present  knowledge  and  technology  of  composite  materials  as 
possible . 


EFFECTIVE  TREATMENT  OF  MEMBRANE  AND  SHEAR  LOCKING  PHENOMENA  IN  CURVED 
SHEAR-DEFORMABLE  AXISYMMETRIC  SHELL  ELEMENT 


A.  Tessler  and  L.  Spiridigliozzi 
U.S.  Army  Materials  Technology  Laboratory 
Watertown,  Massachusetts  02172 

Wide  varieties  of  structural  shells  exhibit  geometric  and  mate¬ 
rial  symmetry.  As  such,  they  are  best  modeled  by  appropriate  shell  of 
revolution  finite  elements  which  are  intrinsically  one-dimensional. 
Thus,  computationally  efficient  analyses  are  possible  owing  to  the 
small  bandwidth  and  dramatically  reduced  number  of  degrees-of -freedom, 
as  compared  to  the  alternative  discretizations  with  more  'expensive' 
two-dimensional  shell  elements. 


Design  requirements  for  thin  and  moderately  thick  shells  and  in 
particular  those  made  of  laminated  organic  composite  materials  neces¬ 
sitate  accurate  modeling  of  the  shell  geometry  and  proper  account  of 
the  effects  due  to  shear  deformation  and,  in  dynamics,  rotary  inertia. 
To  be  effective,  therefore,  an  element  should  be  able  to  accommodate 
these  higher-order  effects  and  closely  approximate  the  meridian  curva¬ 
ture.  Additionally,  nodal  simplicity  is  commonly  desired,  especially 
if  computationally  involved  nonlinear  analyses  are  warranted. 

A  major  difficulty  intrinsic  to  shear-deformable  bending  elements 
is  due  to  a  mathematical  deficiency  called  ' shear,  locking' ;  (the  term 
refers  to  nearly  trivial  kinematic  patterns  manifested  in  the  thin 
regime  due  to  an  overly  stiff  model).  A  mathematically  analogous 
phenomenon  known  as  'membrane  locking'  is  a  common  hindrance  in  curved 
bending  elements.  Thus,  formulating  a  curved  shear-deformable  element 
requires  simultaneous  resolution  of  the  two  locking  phenomena. 


In  this  paper  we  present  a  simple,  yet  extremely  effective, 
two-node,  displacement-type,  axisymmetric  shell  element  derived  from 
Naghdi's  refined  shell  theory  [1]  with  the  aid  of  the  Marguerre  shal¬ 
lowness  equations  [2].  This  work  is  an  extension  of  the  conical  shell 
element  discussed  in  [3]  and  the  Timoshenko-Marguerre  beam  [4],  Here, 
we  adopt  the  anisoparametric  kinematic  approximations  which  employ 
distinctly  different  degree  polynomials  for  the  kinematic  variables, 
rather  than  uniform  interpolations  as  in  the  conventional 
isoparametric  procedure.  Owing  to  these  approximations,  both  the 
membrane  and  shear  strains,  which  serve  as  penalty  constraints  in  the 
thin-element  regime,  possess  physically  admissible  penalty  modes.  The 
formulation  is  further  enhanced  by  an  innovative  treatment  of  the  two 
penalty  parameters  (nondimensional  multipliers  of  the  shear  and  mem¬ 
brane  strain  energies),  giving  rise  to  the  appropriate  relaxation  of 
the  penalty  straining  for  coarsely  discretized  models.  To  focus  on 
the  core  issues,  only  the  axisymmetric  response  of  linearly  elastic 
isotropic  shells  is  considered;  note  that  extensions  to  asymmetric 
response  and  laminated  composites  are  straightforward  once  the  'lock¬ 
ing  issues'  have  been  conceptually  resolved. 
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The  main  features  of  this  two-node,  six  degrees-of-freedom  shell 

V 

element  are  summarized  as  follows  (refer  to  Fig.  1  for  the  shell 
notation) : 

1.  The  meridian  variation  of  the  element  geometric  description 

is  taken  as  cubic  and  is  given  in  terms  of  two  nodal  slopes  0O  and  px. 

.4? 

which  are  assumed  to  be  small  (i.e.,  the  geometry  is  shallow  in  the 

irf 

Marguerre  sense). 

2.  The  kinematic  variables  u  (membrane  displacement),  w  (trans- 

* 

t  t  «r 

verse  displacement)  and  0  (cross-sectional  rotation  about  the  normal 
to  the  meridian)  are  interpolated  in  the  meridian  direction  by  quar- 
tic,  quadratic  and  linear  polynomials  ,  respectively.  The  quartic  u 

• 

is  coupled  in  terms  of  u^,  w^  and  0^  (i=0,l)  degrees-of-freedom;  the 

quadratic  w  is  coupled  in  terms  of  w^  and  0^  degrees-of-freedom;  and 
the  linear  0  is  solely  dependent  upon  its  nodal  values  0..  In  the 
particular  case  of  a  straight  meridian,  u  reduces  to  an  independent 

A! 

linear  function  in  terms  of  u^  degrees-of-freedom. 
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3.  Because  of  the  specific  choice  of  the  kinematic  approximat- 
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ions,  the  membrane  strain,  bending  curvature  and  shear  strain  acting 
in  the  meridian  direction  are  constant.  Thus,  only  one-point  Gaussian 

quadrature  is  needed  to  integrate  the  corresponding  strain  energies 

exactly. 

4.  Element  stiffness  and  consistent  mass  matrices,  and  consis- 
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tent  load  vector  are  derived  via  Hamilton's  variational  principle 
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where  full  Guassian  quadratures  are  employed  on  all  energy  terms. 

5.  The  shear  and  membrane  relaxation  parameters  depend  exclu- 

$ 

sively  upon  element  geometric  and  material  properties  and  yield  a 
well-conditioned  stiffness  matrix  for  any  shell  thickness.  Important¬ 
ly,  these  relaxation  measures  give  rise  to  notable  accuracy  enhance- 
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ments  in  coarsely  discretized  models. 
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Numerical  studies  are  carried  out  focusing  on  the  element  model- 
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ing  capabilities  in  free  vibration  problems. 
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Fig.  1.  Shallowly  curved  axisymmetric  shell  element. 
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EXTENDED  ABSTRACT 


Rubber  track  pads  used  on  Army  tanks  are  subjected  to  large  loads. 
The  loading  environment  consists  of  repeated  compression  loads  applied  by 
the  road  wheels,  penetration  by  rocks  or  metal  debris,  scuffing  and 
cutting.  Failure  modes  include  chunking,  subsurface  cracking  and  blowout 
from  hysteretic  heating.  It  is  not  practical  to  try  and  model  all  possible 
loading  situations  to  which  a  track  pad  may  be  exposed.  Instead,  to  screen 
new  material  candidates,  laboratory  tests  and  analyses  of  samples  can  be 
made.  These  tests  and  analyses  should  subject  the  material  to  the  various 
modes  of  loading  expected  so  that  materials  which  perform  better  than  the 
current  standard  can  be  identified. 

During  the  past  two  years  MTL  has  started  a  research  program  to 
develop  finite  element  analysis  capabilities  for  determining  the  large 
deformations  and  stresses  in  elastomer  (rubber)  materials  used  in  vehicles. 
The  finite  element  programs  developed  can  be  used  to  analyze  and  design  the 
laboratory  tests  mentioned  above.  In  this  effort,  a  total  Lagrangian 
potential  energy  formulation  with  a  Valanis-Landet  form  for  an  energy 
density  function  was  used  to  model  elastomers.  A  nonlinear  finite  element 
algorithm  was  developed  which  used  bilinear  three  node  triangular  elements 
to  minimize  the  potential  energy.  Elastomer  cylinders  loaded  jui  both  end 
thrust  and  by  a  probe  along  the  cylinder  axis  were  analyzed'"  ’  ’  .  More 
recently,  a  modified  form,  of  the  classical  Mooney-Rivlin  energy  density 
functional  was  derived^  which  allows  for  computation  of  rubber 
deformations  without  the  hydrostatic  pressure  as  an  independent  variable. 
Also,  a  .  quadratic  six  node  finite  element  formulation  has  been 
developed'  .  (2-6) 

The  analyses  of  axisymmetr ic  deformations  of  rubber'  can  be 
modified  to  allow  for  the  computation  of  plane  strain  rubber  deformations. 
In  this  effort  we  show  the  plane  strain  formulation  and  present  numerical 
data  obtained  for  a  two  dimensional  model  of  an  Army  T-156  track  pad.  The 
notation  used  to  describe  the  stretching  of  the  material  is  shown  in  Figure 
1.  Both  three  node  bilinear  and  six  node  biquadratic  elements  are  shown  in 
Figure  2.  One  half  of  a  finite  element  model  for  a  T-156  track  pad  is 
shown  in  Figure  3.  The  energy  density  functions  used  are  as  follows: 

VALANIS-LANDEL 
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w  =  2y  E  A.(lnA.-l)  +  A[ln(A  A  A  )]" 
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MODIFIED  MOONEY-RIVLIN 


w  =  Cl(I-3K1/3)  +  c2(J-3K2/3)  +  |xin2K 


where 


and 


w  =  energy  per  unit  undeformed  body, 

A.  =  the  i'th  stretch  ratio, 

1 

2  2  2 
i  =  a^+a2+a2, 

2  2  2  2  2  2 
J  =  a^a2+a2a2+a2a^, 

2  2  2 
K  =  A^A*. 


p,  A  =  the  Lame  constants, 
material  coeffici* 
energy  density  functional. 


Cj,  c^  =  material  coefficients  from  classical  Mooney-Rivlin 
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Figure  1.  Description  of  material  line  deformation  and  stretch  ratios. 
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Figure  2.  Bilinear  and  biquadratic  elements  used  to  discretize  the 
potential  energy. 


Figure  3.  Plane  strain  T-156  tank  track  pad  finite  element  model. 


ANALYSIS  AND  TESTING  OF  PREDAMAGED  SHELL  STRUCTURES 
Richard  E.  Keeffe 
Sr.  Research  Scientist 
Kaman  Sciences  Corp. 

Colo  Spgs,  Colo  80907 

Eric  Reece 

Member  Technical  Staff 
Sandia  National  Laboratories 
Albuquerque,  New  Mexico  87185 

EXTENDED  ABSTRACT 


A  predamaged  shell  structure  may  undergo  significant  reduction  in  load 
carrying  capacity,  dependent  upon, the  location,  size,  and  shape  of  the  damaged 
area.  Kaman  Sciences  Corporation^  (KSC)  and  Sandia  National  Labs  (SNlA) 
under  direction  from  the  U.S.  Army  Strategic  Defense  Command  are  investigating 
analytic  and  experimental  methods  for  evaluation  of  these  effects.  This  paper 
reviews  some  results  of  work  in  progress  and  planned  structural  test  efforts 
aimed  at  supplying  data  for  analytic  correlations. 

Of  primary  concern  are  thin  walled 
cylindrical  and  conical  shell  frustra  which 
have  locally  damaged  areas  as  indicated  in 
Figure  1.  Note  that  these  damages  are 
idealized  and  that  in  general,  damages  may 
be  irregular  in  nature.  As  shown,  three 
basic  loading  conditions  are  of  interest. 

1.  Axial  Tension 

2.  Axial  Compresison 

3.  Lateral  Bending 

Note  that  two  bending  configurations  will 
be  studied;  one  with  the  damaged  area  in 
tension,  the  second  with  the  damaged  area 
in  compression. 

The  basic  analysis  tool  currently  being  considered  for  these  studies 
involves  application  of  the  nonlinear  finite  element  code  ABAQUS,  developed 
and  maintained  by  HIBBITT,  KARLSSON  and  SORENSON,  INC.  The  code  capabilities 
include  static  and  dynamic  modeling  of  solid,  beam,  and  shell  elements,  and 
provides  accurate  modeling  for  arbitrary  magnitudes  of  displacement,  rotation, 
and  strain.  Preliminary  assessment  of  the  applicability  of  the  ABAQUS  formula¬ 
tion  has  been  applied  to  tension  ancLcompression  test  results  obtained  on 
6061-T6  seamless  aluminum  cylinders'^'.  Slot  damages  were  machined  in  one 


FIGURE  1 

DAMAGED  SHELL  GEOMETRY 


H")  This  work  is  being  performed  under  Contract  DASG60-85-C-0051  for  the 
USASDC  Huntsville,  Alabama.  Mr.  Owen  Bender  (USASDC)  is  the  contract 
monitor. 

(2)  This  Testing  was  performed  at  Southwest  Research  Institute,  San  Antonio, 
Texas  under  subcontract  to  Kaman  Sciences  Corporation. 


41 


side  of  these  cylinders  and  compression  and 
tension  tests  (with  clamped  ends)  were 
performed  to  large  deformation  and  frac¬ 
ture.  Axial  load  vs  total  axial  deflection 
predicted  using  ABAQUS  are  compared  to  test 
results  in  Figure  2.  The  comparison  shown 
was  obtained  using  an  initial  radial 
imperfection  of  .01  inches  (.254  mm)  at  the 
shell  midplane.  Additional  variations  in 
initial  imperfections  are  being  analyzed. 

The  final  deformed  mesh  geometry  (X2.5) 
predicted  by  ABACUS  is  presented  in  Figure 
3(a),  along  with  a  photograph  of  the  post 
test  shell  specimen  (Figure  3b).  In 
general  the  comparisons  between  test  and 
theory  are  promising  for  this  single  shell 
geometry. 

FIGURE  2 

ANALYSIS  TEST  COMPARISON 


FIGURE  3 

SHELL  DEFORMATION 


A  more  ambitious  combined  analysis  and  test  program  involving  two  layered 
conical  frustra  is  currently  underway.  Table  I  presents  a  conical  frustra 
test  matrix  planned  for  FY  1986-87.  The  initial  tests  will  incorporate 
single,  circular  hole  damage  patterns.  Additional  tests  planned  for  1987  will 
include  multiple  hole  patterns  and  non  circular  damages.  Actual  testing  will 
be  performed  by  the  Sandia  National  Laboratories  in  Albuquerque,  New  Mexico. 


TABLE  I 

CONICAL  FRUSTRA  TEST  MATRIX 
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CORROSION  CONTROL  BY  COMPOSITE  MATERIAL  SUBSTITUTION  WITH 
APPLICATION  TO  THE  M-939  5  TON  TRUCK 


PIN  TONG,  Chief  Structures  and  Dynamics  Division 
TRANSPORTATION  SYSTEMS  CENTER 
CAMBRIDGE,  MA  02142 

In  an  effort  to  improve  the  corrosion  resistance  of  the  5-ton  M-939  truck  and  to 
lighten  its  weight,  the  Army  is  evaluating  the  replacement  of  structural 
components  which  have  high  potential  for  environmental  deterioration  with 
equivalent  composite  components.  Asa  part  of  this  overall  effort,  several  graphite 
epoxy  chassis  were  fabricated  by  Ewald  Associates.  A  schematic  of  the  chassis  is 
shown  in  figure  1.  These  chassis  were  designed  to  match  the  strength  and  stiffness 
characterisitcs  of  the  steel  chassis  and  weigh  about  50%  less.  One  of  the  chassis 
has  been  installed  in  an  actual  truck  and  is  being  used  for  demonstration  purposes. 

The  present  study  involves  the  evaluation  of  the  prototype  chassis  both  as 
fabricated  and  as  installed  in  the  5-ton  truck  for  corrosion  resistance,  structural 
integrity,  and  damage  tolerance.  Of  interest  is  the  sensitivity  of  the  composite 
frame  mechanical  details  to  stress  concentrations,  fracture,  shock  and  vibration 
and  environmently  produced  degradation  particularly  from  moisture  and  abrasion. 
Additionally,  other  structural  components  which  are  critically  sensitive  to 
corrosion  and  may  afford  a  great  weight  saving  potential  will  be  identified  and 
alternative  composite  designs  developed.  Finally,  a  prototype  truck  employing  as 
many  light  weight  corrosion  resistant  components  as  possible  will  be  designed, 
fabricated,  and  tested. 

In  the  initial  phase  of  this  project  a  methodology  has  been  developed  to  evaluate 
the  composite  frames.  This  methodology  involves  identifying  stress  sensitive 
locations,  characterizing  static  strength  and  dynamic  flexibility  of  the  composite 
frame,  and  evaluating  the  in  service  stresses.  The  first  step  has  been  to  compare 
the  stiffness  characteristics,  as  well  as  the  first  bending  and  warping  natural 
frequencies  of  the  composite  frame  with  the  steel  frame.  In  addition,  the  worst 
case  loading  conditions  have  been  identified  and  a  simple  dynamic  model  proposed. 
Estimates  are  developed  for  the  effect  of  the  assembly  modifications  made  to  the 
Ewald  frame  currently  installed  in  an  actual  truck. 

45 


Test  plans  for  the  truck  with  the  composite  chassis  include  estimating  the  principle 
connection  stiffness  and  damping  characteristics  of  an  assembled  truck.  These 
parameters  are  required  for  the  dynamic  model  of  the  truck  in  conjuction  with  road 
spectrum  data,  for  the  prediction  of  life  cycle  estimates.  Of  critical  concern  is 
the  potential  for  environmental  deterioration  of  the  frame  especially  at  or  near  the 
mechanical  fasteners.  Tests  conducted  on  coupons  or  sample  mechanical  details 
are  planned  to  evaluate  effects  such  as  moisture  content  and  operating 
temperature  (-25°F  to  +  125°F).  These  tests  will  include  sensitivity  to  freeze- thaw 
cycles  and  mechanical  abrasion. 

Numerous  truck  components  have  been  identified  as  corrosion  sensitive  with  some 
potential  for  weight  reduction,  also,  designs  for  these  components  will  be 
developed  and  evaluated  for  structural  capacity,  particularly  in  reference  to 
abrasion,  environmental  deterioration,  and  damage  tolerance.  The  final  phase  of 
the  project  will  be  the  fabrication,  assembly  and  service  testing  of  a  prototype 
truck  encompassing  all  the  light  weight  corrosion  resistant  components  developed. 


Figure  1.  CURRENT  STEEL  5-TON  TRUCK  CHASSIS 
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Many  high  strength,  low  density  powdered  metals  are  being 
developed  for  aerospace  applications.  Advantages  in  mechanical 
properties,  such  as  lower  density,  higher  fracture,  toughness,  and  much 
higher  strength,  in  powdered  metallurgy  wrought  products  are  due  to 
the  fine  grain  structure  resulting  from  rapid  solidification  achieved 
by  powder  atomization.  Powdered  metallurgy  (P/M)  offers  signficant 

improvements  in  material  properties  over  ingot  metallurgy  (I/M). 

/ 

Although  the  physical  and  mechanical  properties  of  P/M  materials 
are  still  being  quantified,  properties  of  P/M  alloys  are  promising 
greater  load-carrying  capability  and/or  reduced  weight  in  high-stress, 
fracture  critical  airframe  applications.  Even  with  this  promising 
outlook,  the  applicabilty  of  powdered  metals  in  fracture  critical 
aerospace  components  is  uncertain.  In  order  to  help  maximize  the 
performance  of  military  aircraft,  this  research  program  is  being 
conducted  to  investigate  the  crack  growth  characteristics  of  a 
representative  powdered  aluminum,  and  to  characterize  the  initial 
quality  and  fatigue  crack  growth  rates  of  P/M  aluminum  for  durability 
and  damage  tolerance  analysis.  In  addition  to  determining  if 
traditional  crack  growth  fracture  mechanics  analysis  is  possible,  this 
program  will  ensure  that  structures  made  of  P/M  aluminum  can  be 
designed  to  ensure  structural  integrity. 

This  research  program  investigates  the  characteristic  properties 
of  extruded  bar  stock  7091  T7E69,  a  first  generation,  high  strength 
P/M  aluminum  produced  by  ALCOA.  The  test  specimens  of  dimensions 
M " x . 5 " x 4 8 "  have  several  different  flaw  geometries  including  elox 
center  starter  notches  of  approximately  .1"  (figure  1A),  open  holes 
0.25"  in  diameter  (figure  IB),  open  holes  0.25"  in  diameter  with 
corner  flaws  (figure  1C),  and  open  holes  0.25"  in  diameter  with  thru- 
the-thickness  flaws  (figure  ID).  The  flaw  geometries  were  prepared 
sixteen  inches  (16")  apart  to  ensure  no  load  interaction  between  the 
two  holes  or  notches. 


The  test  plan  includes  four  tensile  tests  to  determine  the 
mechanical  properties  of  7091,  and  ninety-four  fatigue  crack  growth 
tests  of  constant  amplitude  and  variable  amplitude  load  histories. 
The  variable  amplitude  load  history  being  used  is  a  programmed  F-16 
load  history  containing  predominately  tension-tension  cycles  as  well 
as  a  few  compression  cycles  with  each  block  representing  400  flight 
hours . 

To  date,  all  tensile  tests  and  twenty-five  fatigue  tests  have  been 
completed.  The  fatigue  test  specimens  show  extreme  out  of  plane 
cracking  and  tunneling  (figure  2).  Believing  that  the  out  of  plane 
crack  growth  was  due  to  the  thickness  of  the  specimen  (.5")  and 
residual  surface  stresses,  the  specimens  were  modified.  In  an  attempt 
to  alleviate  surface  stresses  and  to  control  the  non-planar  crack 
growth,  the  specimens  were  milled  to  3.9"x.4"x48".  Then  the  maximum 
stress  level  was  reduced  from  30  KSI  to  25  KSI.  In  the  thinner 
specimens,  the  crack  growth  rate  decreases,  and  the  total  life  is 
increased  by  twenty  six  blocks  on  the  average.  Also  the  angle,  G,  at 
which  the  crack  deviates  from  the  plane  of  the  flaw  is  decreased  by 
forty  percent  on  the  average  (figure  3).  The  thinner  specimens 
appear  to  exhibit  better  crack  growth,  yet  this  crack  growth  still 
exhibits  extreme  non-planar  growth  and  tunneling. 

The  fatigue  crack  growth  of  P/M  7091  is  extremely  consistent.  The 
characteristic  nature  of  non-planar  crack  growth  and  tunneling  is 
present  in  all  of  the  fatigue  test  specimens.  The  non-planar  crack 
growth  and  tunneling  is  not  dependent  on  the  thickness  of  the 
specimen  or  the  maximum  stress  level.  For  this  reason,  it  is  believed 
that  the  non-planar  crack  growth  is  completely  normal  for  P/M  7091  and 
that  linear  elastic  fracture  mechanics  may  not  apply  to  P/M  metals.  A 
new  method  must  be  developed  to  analyze  the  fatigue  crack  growth  of 
P/M  7091  and  all  other  P/M  metals. 

Tensile  tests  showed  that  P/M  7091  demonstrates  an  average  yield 
strength  of  72.7  KSI,  an  average  ultimate  strength  of  83.9  KSI,  and  an 
average  modulus  of  elasticity  of  9-9x10°  PSI.  I/M  7075  T7  demonstrates 
a  yield  strength  of  59  KSI,  an  ultimate  strength  of  80  KSI,  and  a 
modulus  of  elasticity  of  10.4x10°  PSI.  Therefore  P/M  7091  has  an 
improved  yield  strength  and  ultimate  strength  over  I/M  7075,  but  a 
decreased  modulus  of  elasticity. 

Before  any  more  fatigue  testing  is  performed,  life  predictions 
will  be  developed  for  P/M  7091  using  the  Walker  equation.  Then 
testing  will  be  resumed  to  check  the  accuracy  of  the  developed 
equation.  Adjustments  to  the  equation  will  be  made  if  necessary. 

Conclusions  of  this  research  on  the  durability  of  powdered 
aluminum  structures  will  be  documented  when  the  fatigue  test  program 
is  completed. 
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The  problem  of  erosion  and  wear  of  gun  barrels  has  long  been 
recognized.  Ceramic  materials  offer  the  most  favorable 
combination  of  high  hardness,  compressive  strength,  erosion 
resistance  and  chemical  inertness  at  high  temperature.  In 
recognition  of  the  potential  of  ceramics  for  gun  barrel  life 
improvement,  MTL  and  DARPA  contracted  with  SACO  Defense  Systems 
Division  to  investigate  the  feasibility  of  a  ceramic-lined 
barre  1 . 

This  work-in-progress  report  describes  the  analysis  and 
failure  evaluations  of  the  tests  performed  on  three  potential 
ceramric  liner  materials.  The  ceramic  lined  gun  barrel  test 
assembly  is  shown  in  Figure  1.  The  ceramic  liners  have  an 
interference  fit.  They  are  assembled  by  a  shrink  fitting  process 
first  to  the  steel  sleeve,  then  to  the  steel  jacket.  Shrink 
fitting  is  performed  in  order  to  introduce  favorable  compressive 
stress  in  the  ceramic,  which  could  negate  any  tensile  stresses 
arising  during  firing.  Finally,  the  liner  retainer  and  nozzle 
attachments  are  assembled.  Following  that  the  ceramic  liner  is 
thoroughly  inspected  using  a  horoscope  before  the  firing  tests. 
After  each  firing  the  ceramic  liner  is  also  inspected  for  any 
cracks . 

The  three  ceramic  materials  tested  were«i-siC,  Si^N^  and  PSZ/TS. 
Table  1  gives  a  summary  of  the  firing  test  results  and  failure 
modes  observed. 

In  order  to  understand  the  cause  of  failure  in  each  case,  a 
finite  element  analysis  was  performed.  The  finite  element 
analysis  attempted  to  simulate  the  shrink  fit  operations  and 
firing  conditions.  It  should  be  mentioned  that  SACO  performed  a 
Lame  analysis  in  order  to  evaluate  the  interference  fit  and 
firing  stresses,.  The  finite  element  model  in  Figure  2  emphasized 
the  details  at  muzzle  end  of  the  ceramic  liner  where  most  of  the 
failures  occurred. 

Figures  3  to  8  show  the  hoop  and  axial  stresses  due  to  shrink  fit 
and  firing  and  are  used  in  the  failure  analysis. 


*This  program  is  Co-sponsored  by  MTL  and  DARPA. 


a 


.'i. 


s 


3 

I 


The  presentation  will  discuss  the  analysis  results  for  the 
assembly  operation  and  firing.  The  failures  observed  of  the 
three  materials  are  correlated  with  the  analysis  results.  It 
will  also  cover  the  use  of  the  finite  element  results  in  a 
Weibull  analysis  to  evaluate  the  probability  of  failures. 

Dynamic  effects,  design  changes  and  recommendations  for  improved 
performance  will  be  discussed  also. 
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FIGURE  2  -  Finite  element  model  of  gun  barrel  section 
containing  ceramic  liner 


FIGURE  3  -  Contour  plot  of  hoop  stress  FIGURE  4  -  Contour  plot  of  axial  stress 

due  to  shrink  fit  and  firing  due  to  shrink  fit  and  firing 

with  ano<-SiC  liner  with  an  °^-5iC  liner 
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